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SUMMARY 

A statistical theory based on the presence of small crystalli- 
zation nuclei suspended in water is developed to explain experimental 
results showing that on the .average small droplets can he supercooled 
to lower temperatures than large ones. Small nuclei of crystalliza- 
tion are assumed responsible for causing supercooled water to freeze 
spontaneously. 

The average behavior of supercooled droplets is reproduced on the 
basis of probability theory with an assumed distribution of crystalli- 
zation nuclei with respect to the temperatures at which the nuclei 
cause freezing. The most probable distribution curves of spontaneous 
freezing temperatures for water droplets of various sizes within the 
size range found in clouds are obtained. 


INTRODUCTION 

The prolonged operation of unprotected aircraft in supercooled 
clouds has been found to be a hazard. For the icing hazard to be over- 
come, the physical conditions under which supercooled clouds exist must 
be understood. The physics of supercooling and freezing is incompletely 
understood, as shown by conflicting opinions expressed in the litera- 
ture. The following factors have been described as causes of freezing 
in supercooled water samples: (1) introduction of ice crystal; (2) mech- 

anical shock, such as vibration, scratching of solid surfaces together 
in water sample, and sound waves; (3) length of time water has existed 
in supercooled state; (4) impurities; and (5) droplet size. 

Investigators seem to agree that an ice crystal introduced into a 
sample of supercooled water will cause the water to freeze. Recent 
experiments by Dorsey (reference 1) indicate that mechanical shock, 
however, is of little importance in freezing of supercooled water; like- 
wise, time apparently is not a factor. From the extensive investigation 
of reference 1, it is concluded that small solid particles suspended in 
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the water or irregularities on the container surface . cause the super- 
cooled water to freeze at a given temperature . These hypothetical par- 
ticles are called motes and the temperature at which a sample freezes 
without any apparent, external stimuli is called the spontaneous freezing 
temperature. The spontaneous freezing temperature was found in most 
cases to he a characteristic of each sample, which led to the mote 
hypothesis . 

The observations described in reference .1 have been substantiated 
somewhat in reference 2 in which a statistical investigation of water- 
droplet size as a factor in determining the spontaneous freezing tem- 
peratures of the droplets is presented. By observation of a large number 
of droplets, the following results were obtained: (l) The frequency- 

distribution curve with respect to spontaneous freezing temperature of a 
given droplet size was found to have a marked peak or mode; and (2) the 
average spontaneous freezing temperature of each size group decreased 
with decreasing droplet size. 

In reference 3, the conclusions are made from observation of a 
number of droplets that the spontaneous freezing temperature of droplets 
from 400 microns down to approximately 50 microns rapidly decreased with 
size and that droplets larger than 400 microns had a practically constant 
freezing temperature. The results of reference 3 are qualitatively 
verified in reference 2 but only to the extent that the average freezing 
temperature for a given size was found to be a function of the size. 

Similar results of a qualitative nature have been obtained with 
substances other than water. In reference 4, the possible degree of 
supercooling in sulphur droplets was found to increase with decreasing 
droplet diameter. In addition, it was found more difficult to super- 
cool sulphur droplets on a metal surface, such as brass or aluminum, 
than on a glass surface. -The effect of droplet size on spontaneous 
freezing temperature was ascribed to the fact that smaller droplets are 
less susceptible to external influences than larger droplets. The sur- 
face was thought to have an influence on the crystallizing forces within 
the supercooled droplet. Substances such as gold, platinum, rhodium, 
mercury, and other metals have been supercooled to a greater extent in 
the form of small particles than in bulk (reference 5) . The effect of 
droplet size on spontaneous freezing temperature is ascribed in refer- 
ence 5 to the presence of small particles of impurity in the liquid. 

When the liquid is broken up into small droplets the foreign particles' 
are isolated in a few drops allowing the remainder to supercool to the 
temperature at which molecules of the liquid combine to form crystals 
without the aid of motes. 
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The hypothesis of reference 5 comhined with the mote concept of 
reference 1 indicates a way of explaining the statistical results pre- 
sented in reference 2. A statistical theory that accounts for the two 
conclusions of reference 2 on the basis of motes or any other agent 
causing nuclei of crystal growth to form was developed at the NACA Lewis 
laboratory and is presented herein. 


THEORY 

An assumption is made that a large number of motes are present in 
liquid water. Also, each mote is assumed to be associated with a defi- 
nite spontaneous freezing temperature. The freezing temperature of a 
water sample is governed by the mote in the sample that is associated 
with the highest freezing temperature. 

In the experiment reported in reference 2, water was condensed 
from the atmosphere in the form of frost on a polished metal surface 
by lowering its temperature to approximately -35° C. The frost was sub- 
sequently melted and droplets formed on the surface . Repeated applica- 
tion of this procedure made a large number of droplets of various sizes 
available for determining spontaneous freezing temperatures. Because 
the water condensed from the atmosphere was. a very small fraction of 
the total amount of water available in the form of vapor, the process 
was practically the same as drawing small samples of volume AV from 
a large reservoir of liquid water of volume V. Because motes of all 
kinds are assumed to be present everywhere, it is reasonable to con- 
sider them to be present in the same relative proportions in the hypo- 
thetical volume of liquid water V. 

The problem to be solved is the determination of the most probable 
frequency-distribution curve of spontaneous freezing temperatures sub- 
ject to the following conditions: 

(1) There are N 0 , H x , N 2 , K 3 , . . . N T , . . . motes of kinds 

0, 1, 2, 3, . . . T, . . ., respectively, present in the volume V, 
where T is the absolute value of the temperature in °C at which 
the mote of kind T causes water to freeze. 

(2) The small volumes that are drawn from V are of volume AV, 
and the total number of volumes AV obtainable from V is V/AV = p. 

(3) The number of draws in each set of draws from V is so small 
that p is practically unchanged while a set is being drawn. 
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First the motes of kind. T are considered. There are Nrp motes 

available for distribution among the p small volumes AV, into which V 
may he divided. The total number of ways of distributing Mp objects 

among p cells is p 1 ^. Each way of distributing Nip objects may be 
called a complexion. Next, regarding the group of all possible com- 

N +1 

plexions , there are p cells in each complexion and therefore t p 
cells in all the complexions. The probability of drawing a small volume 


in one draw with no motes of kind T in it is , > where Q is the 

Nm+X 

Nir+1 ^ 

number of cells that among the p 1 cells in all the possible complex- 
ions contain no motes. The probability of drawing a small volume contain- 


ing one or more motes is 1 


Q 

Nip+1* 


A simple example should serve to make the foregoing logic clearer. 
Two motes are supposed to exist in a volume of water V. From V a 


volume AV = i Y ' is to be drawn. The probabilities that in a single 

w 

draw one or more motes will be in the half volume drawn or that no motes 
will be contained in it are to be determined. For solution of the prob- 


lem, the volume Y may be schematically represented by two squares of 
equal, area side by side. Thus, if the motes are numbered 1 and 2 


for identification, the complexions may be pictured as follows: 



The number of possible- complexions is = 2^ = 4 and the total number 

of volumes or cells in all the complexions is ^T+l _ ^5 = 8. The num- 
ber of cells containing no motes is two and the number containing one or 
more motes is six. The probability of drawing a cell with no motes in it 


is therefore 



and with one or more motes in it, 1 


Q 

Nrn+1 

P ■ L 


In the preceding example a particular case has been evaluated. 
General expressions for the probabilities must now be obtained in terms, 
of , p and %>. As in the previous example, cells are lined up side by 
side numbered from 1 to p. In a given complexion there are Nip p, 

Nt, 2> • • • ^T,p motes in the cells 1, 2, . . . p, respectively, where 
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the second subscript is the cell number and N T j_ is the number of motes 
in the cell: 


1 

2 

3 

4 

5 

N T,1 

W T,2 

W T,3 

N T,4 

%,5 


(Mr 


TJJ 


m 


SiHl 


The total number of motes causing freezing at a temperature of -T° C is 


'% = W T,1 + N T,2 + 


+ 


.Na 


The total number of complexions having Ny ± equal to 0 is (4-1) , 

because Wj objects are restricted to distribution among 4-1 cells 

N 

by the condition ^=0. Likewise, there are (4-1) T complexions 


subject to each of the conditions 2 =°# 


%,i=°> 

Nm 


W T,4 -0 ' 


respectively. Thus, there are altogether 4(4-1) cells among the 


Krr 


n t +i 


4 * complexions that contain no motes. But there is a total of 4 
cells In all the possible complexions. The probability of drawing a cell 


with no motes in it is therefore 
one with one or more motes in it is 


4 N 


and the probability of drawing 


/ _i|t 

VT7. 


In order to complete the solution of the problem with all kinds of 
mote present, the probability P(T,|a) of drawing a cell with one or more 
motes causing freezing at a given temperature and none causing freezing 
at a higher temperature must be found. The probability of drawing a 
cell that would freeze at a temperature of -T° C is equal to the product 
of the probability of drawing a cell with no particles of kinds 
0, 1, . . . T -1 and the probability of drawing a cell with one or more 
particles of kind T, that is, 


P(T,n) 



Nq+N^+N2+ 



a) 


The most probable frequency -distribution curve of spontaneous freezing 
temperature for a given value of 4 is given by equation (1) . 
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•Because the quantity p = v/AV contains "both AV and V, the 
question arises as to the possible dependence of P(T,u) on V as well 
as AV. It can be. shown that for AV/V small enough so that terms of 
the second or higher order may be neglected, the function P(T,p) is 
independent of V. Samples of volume AV are drawn, from two large 
volumes V and V* such that AV/v and AV/V' are very small. The 
mote densities (numbers per unit volume) of various kinds are the same 
in both volumes, namely Pq, P]_, Pg> • • • . The corresponding expres- 

sions for P(T,|i) are 




( p 0 +p l +p 2 + • • • + Pp_i)V 




and 


P* .= 1 - 


( p 0 +p l +p 2 + * * • + 


- i-ff 


With only the first two terms in the hinominal expansion of 

/ \( p o +p i +p 2 + ? ‘ * + p T )V / \ p T V 


— J , it is found that 


P = P’ 

which proves P(T,p) to be independent of the choice of V provided 
AV/v is small enough. Even the higher-order terms are practically • 
independent of V if the exponents are large compared to 1. (See the 
appendix.) 

With the expression for P(T,p) formulated in equation (l), a mote- 
density function of freezing temperature must be found such that the 
resulting probability function agrees at least in its general character- 
istics with the distribution curves presented in reference 2. The 
"fch 

(T+l) term of the geometric progression 

a, ar, ar 2 , . . . ar^ 1 , ... 

where a is the number of motes in the volume V causing freezing at 
0° C and r is an arbitrary constant, was found to be a suitable mote 
frequency -distribution function. Thus the sum of all motes causing 
freezing between 0° C and -T°.C in V is 
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The following fora of equation (1) is obtained by substitution of 
the foregoing expression: 


p »» - (hr) 

■because 


N 0 +N l +N 2 +N 3 + * • 

The nature of equation (2) is evident in figure 1 where curves corres- 
ponding to particular values of p and Pq and various values of r 

are plotted* The probability curves all have a max i mum and are some- 
wnat unsymmetrical • Furthermore , as r increases the maximum value 
of P increases and the peak of the curve becomes narrower. The curves 
of figure 1 are comparable in form to the frequency-distribution curves 
of reference 2. 

In order to fit the experimental average freezing temperature - 
droplet size curve with the theoretical curve, the mode (or maximum) of 
P(T,p) is found by taking the derivative dP(T,p)/dT and solving the 
equation 

. o 

dr 

for T. The mode or the maximum of the probability curve, which is 
almost symmetrical, is close enough to the average so that they may be 
considered equal (table I). 

By performing the foregoing operations, the following equation for 
the mode % as a function of p, a, and r is obtained: 
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Equation (3) is now reduced to a simple approximate form. Because V 
is much larger than • AV, second- or higher-order terms in AV/v are 
neglected with the result that 


J£_ 

u-1 



1 


+ 


AY 
V ' 


(4) 


and 


log -H— ~ log (l + —) = 0.4343 — 

H-l V V 


(5) 


In the experimental curve of reference 2, is plotted as a function 

of the diameter D of a hemispherical drop of volume AV so that for 
comparison of theory with experiment equation (4) becomes 


V 

H-l 



or 

log S 0.4343 

U-1 

Furthermore, the mote -density function is ^ r T , and the mote- 

8r 

density constant is Pq = — , Thus, by substitution of equation (6) in 
equation (3) 


«D° 

12V 


( 6 ) 


t m = ‘ 


log (log r) - 3 log D + log 
log r 


8.78 

PO 


(7) 


is obtained. As a result of the assumed mote -density function, is 

therefore to the first order of approximation a linear function of log D. 
The. linear relation between T M and log D is a basic one that can readily 

be compared with the experimental data of reference 2 . 
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The slope of the linear equation (7) is 


T 2,M~ T 1,M _ _ 5 

log Dg “ log Dq log r 


( 8 ) 


where the subscripts 1 and 2 indicate two distinct values of T^ corres- 
ponding to two droplet sizes Dq and D 2 , respectively. The constant po 
is found from equation (7), after setting %•= 0> to be 


PO 


8,78 log r 

D 0 3 


(9) 


where Dq is the value of D when Tm = 0. From equation (8), the slope 
of equation (7) is apparently a function of r alone, and pq according 
to equation (9)- determines the intercept. Equations (8) and (9) will be 
subsequently useful in evaluating the constants r and pq ’ from the 
slope and the intercept of the experimental average temperature -log D 
curve . Then with the constants pq and r known', the theoretical 
probability curves may be computed from equation (2) and compared with 
the experimental frequency-distribution curves of . reference 2. 

In addition to equations (8) and (9) defining the mode and the inter- 
cept, an approximate transformation equation is derived in the appendix. 
With this transformation, a probability curve for one set of values of 
the parameters pq and AT can be transformed simply to a curve corres- 
ponding to another set of values. If r is kept constant and it is 
desired to transform from a distribution curve for a volume AVq and 
mote-density constant Po,i to another for AV2 and Po,2> "then for a 
given value of the probability P the value of the freezing temperature 
Tg corresponding to AVg and Po,2 ts rented to Tq corresponding 

to AVq and Po,i by 

log (p 0 iAV 1 /po 2 AV 2 ) 

T 2 = Tq + 2 1 (10) 

log r 

This transformation is possible, because the family of curves correspond- 
ing to all values of AV and Pq are identical in size and shape but 
are displaced in the T coordinate. (See the appendix.) 
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A method of transformation from a curve for rj_ to another for rg, 
keeping Pq and AY constant, is unavailable, but a set of curves 
for various values of r is given in figure 1 from which curves for 
other values of AV and pq may be computed using the preceding trans- 
formations . 

RESULTS AND DISCUSSION 

/ 

For comparison of the theory of this, report with the experimental 
results of reference 2, the data of reference 2 for the average freezing 
temperature -size curve are replotted on semilog paper, and an irregular 
curve that can be fitted reasonably well with a straight line is obtained 
(fig. 2). This result agrees with the linear relation between Tm and 
log D indicated by equation (7). The slope of the. straight line, as 
drawn in figure 2, is -7.21 and the value of Dq is 24 (cm). From the 
slope, the value r is found by equation (8) to be 2.60 and, from the 
intercept Do and equation (9), Pq has the value 2.64x10-4 (cm-3) . 

The linear equation that best fits the data of reference 2 is therefore 

,Tj 4 = 9.95-7.21 log D 

With the constants r and Pq now known, it is possible to com- 
pute the probability curves for various droplet sizes. A comparison 
between these curves and the corresponding experimental frequency- 
distribution curves of reference 2 is given in figure 3 where histograms 
of the experimental data are reduced to the probability scale . Fig- 
ures 3(b) to 3(d) represent the results obtained for the dx.oplet sizes 
with the greatest amount of experimental data, and figures 3(a) and 3(e) 
represent the results for the extremes of the size range investigated 
where the data were limited. In the computation of the ‘probability 
curves of figure 3, a convenient value of V = 10 5 (cm 3 ) was assumed. 

The value of V does not influence the curves if it is large compared 
with AV. 


The agreement between the theoretical and experimental curves in 
figure 3 is qualitatively very good, but the following differences in 
detail exist. The modes of the theoretical curves differ from those of 
the experimental curves corresponding to the deviation of experimental 
points from the straight line drawn through them in figure 2 . The 
experimental probability curves are broader and have less pronounced 
peaks than the theoretical curves. The standard deviations (root mean 
squares of the deviations from the mean freezing temperature ) of the 
theoretical curves are constant and equal to 1.38° C, but the standard 
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deviations of the experimental curves are between 2 and 4 times as great 
and are roughly oonstant in the size range below 100 microns (fig. 4). 
Above 100 microns there is a sharp increase in the standard deviation of 
the experimental data to a higher approximately constant value of 3.8. 

The foregoing difference between the experimental and theoretical 
curves cannot be ascribed to the fact that the experimental data are 
grouped into narrow droplet-size ranges to provide the freezing- 
temperature distribution curves for a given average droplet size. A 
composite theoretical curve formed by combining theoretical curves from 
the middle and the two extremes of the experimental -size group differed 
very little from the unmodified theoretical curve, as shown in fig- 
ure 3(c) . 

Better agreement between theory and experiment is obtained if 
several different mote -density functions are postulated and the corres- 
ponding theoretical curves combined. This procedure may be justified, 
because the experimental data were taken on several days distributed 
over a period of three months. Data taken each day covered almost the 
entire range of droplet sizes investigated. If impurities were the 
active agent, the mote -density function could change from day to day 
during that time. Thus, the combination of data taken at different 
times could easily broaden the experimental distribution curves as com- 
pared to the theoretical curves . Changes in r and pq could modify 
the slope and the intercept, respectively, of the % - log D curve 
(fig. 2). 

The modes of the theoretical distribution curves may be shifted 
without changing the slope of the Tj*j -log D curve by varying the 
constant Pq. The weighted combination of several curves corresponding 
to values of pq between 10"5 and 10“3 (cm"3) (see table II) results 
in theoretical curves that fit the experimental curves much better than 
the unmodified theory, as shown in figure 3. The modified theoretical 
curves, however, have a constant standard deviation equal to 2.12° C, 
which still does not agree with the standard deviations of the experi- 
mental curves, as shown in figure 4. Better agreement might have been • 
obtained by varying the values of r as well as of Pq; but the 
additional work involved in varying r was not considered worthwhile, 
because the variation of Pq was only postulated to demonstrate the 
plausibility of the hypothesis concerning the behavior of the experi- 
mental standard deviations. The good agreement obtained between theory 
and experiment indicates that .variation of the mote -density function 
among several groups of combined data can account for the discrepancy 
between theory and experiment. 
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Deviations from the most probable curves may by chance occur in the 
observed distribution curves, but deviations of a given magnitude become 
less probable as the number of observations increases . This principle 
is illustrated by the greater regularity of distribution curves based on 
a large number of observed droplets as compared with those based on a small 
number of droplets. For example, the curves in figures 3(a) and 3(e) 
for which 39 and 35 droplets, respectively, were observed are quite 
irregular as compared with the curve in figure 3(c), for which 692 drop- 
lets were observed. 


CONCLUDING REMARKS 

In the derivation of the foregoing statistical theory, emphasis' has 
been placed on the mote hypothesis of Dorsey, but the theory is' also 
applicable to homogeneous nucleation as described by Volmer (reference 6). 
According to reference 6, the random small-scale fluctuations in an 
impurity-free supercooled liquid may cause crystal nuclei to form, result- 
ing in crystal growth from the mother liquid. The presence of motes in 
the liquid would mask the homogeneous nucleation. If it were possible 
to remove all foreign nuclei, then the mote distribution function Nq> 
would merely be renamed the homogeneous nucleus distribution function 
and the logic of the theory in this report would apply as in the case 
of the motes except that Nj would now be a. function increasing with 
time. Thus, a container of liquid kept at a constant temperature of 
supercooling should freeze after a period of time. But the experience 
that Dorsey (reference .1) had in maintaining three test tubes' at a tem- 
perature between -8.0° and -10.3° C for 312 days might seem to invalidate 
the homogeneous -nucleation theory. The apparent inconsistency that 
exists between the results of reference 1 and those expected from the 
homogeneous -nucleation theory may possibly be resolved by the following 
logic. The temperature range in which the experiment of reference 1 
was performed may be one of low nucleus -formation- probability even over 
a period of 312 days. But if the probability of nucleus formation 
increases very rapidly at some temperature much lower than -10° C, the 
most probable length of time that a given water sample can remain at a 
low constant temperature of supercooling may be sufficiently reduced to 
make maintenance of the water sample in the supercooled state for a very 
long time practically impossible. 

If the spontaneous freezing temperatures of droplets are princi- 
pally determined by the presence of motes, the results of this report 
can be extended to supercooled clouds, because there should be no 
essential difference between the droplets supported on a surface and the 
droplets suspended in the atmosphere except for possible variations in 
the mote distribution curves as a function of time and place. The 
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results of this report indicate that no well-defined temperature exists 
at which supercooled cloud droplets are transformed into ice particles. 
As the temperature is decreased in a given cloud, more and more droplets 
will tend to freeze until a point is reached at which very few droplets 
are left unfrozen. For example, in a cloud of uniform droplet size 
with hemispherical diameter D of 34.5 microns or an actual spherical 
droplet diameter of 27.4 microns which is fairly common in clouds, prac- 
tically all the droplets would he frozen at a temperature of -34° C . 

On the other hand, practically all the droplets in the cloud would 
remain unfrozen down to a temperature of -18° C and only about half 
would he frozen at a temperature of -28° C. In clouds with a droplet 
size of 5 microns or greater, therefore, very few supercooled droplets 
exist at temperatures lower than -35° C and the existence of an icing 
cloud at such low temperatures is . prohahly an unusual occurrence . 

The relative proportions of frozen and supercooled droplets in a 
cloud are important in determining the stability of the cloud. A cloud 
is considered to he stable if factors causing rapid dissipation of the 
cloud are absent. One of the factors that may cause a supercooled 
cloud to become unstable is the existence of ice particles in the cloud. 
Ice particles in a supercooled cloud have a lower vapor pressure than 
the surrounding supercooled drops. As a result of the vapor -pres sure 
gradient existing between the liquid droplets and ice particles, water 
tends to evaporate from the droplets and condense on the ice particles . 
Thus, if ice particles are present insufficient number distributed 
evenly throughout the cloud, the ice particles will grow at the expense 
of the. supercooled drops until they fall cut of the cloud in the form 
of precipitation (reference 7). If too few ice particles appear in a 
supercooled cloud, they will grow and fall out of the cloud leaving it 
practically intact. On the other hand, if too many ice particles appear, 
the cloud is merely transformed to an ice -particle cloud. The presence 
of ice crystals in supercooled clouds is considered to be the principal 
factor initiating precipitation. 

Determination of whether the statistics of freezing are affected 
by change of environment requires further investigation. For example, 
if impurities are the cause of freezing in supercooled water, then the 
average freezing temperature of a given droplet size should be depressed 
by maintaining as high a degree of cleanliness as possible . In addition, 
if cleanliness is found to have the anticipated influence, the effect 
of adding finely divided materials of kpown constitution to water can be 
investigated. In the experiments on supercooled water that have been 
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conducted up to the present time, the water has always "been in contact 
with a container or support. In order to attain conditions that approach 
those in clpuds, a method of investigating the free zing- temperature 
statistics of unsupported droplets must he found. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, May 3, 1950. 
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APPENDIX - DERIVATION OF TRANSFORMATION EQUATIONS 

Equation (10) for transformation of a probability curve 
corresponding to one set of values of the parameters AV and p 0 to 
another set of values is derived from equation (2) written in the form 




P(T,AV,P 0 ,r) = (l “ f) 0 r_1 1 “ (* - f) 


P0^ v 


( 12 ) 


Equation (12) must first he expanded in a series to obtain an approxi- 
mate expression, which is used to show the dependence of P on AY and 

' *vVo(# 


Pq. The members of equation (12), (l - 

\ v 

are written in a binomial expansion as follows: 


and 1 - 


AV 

TV 


Vo rTv 


AV) nV , _. Tr nV(nY-l) AV 2 nV(nV-l)(nV-2) AY 3 
V/ - lrnAV + 2 ! v 2 3 r " v 3 + 


where n represents 




or PqtT. Because V is very large, 


nV in general is so large that nV-1, nV-2, . . . are practically equal 
to nV so that the binomial expansion may be written as 


(l - Yj Ul ~ 1-nAV + AV 2 - —■ AV 3 + . . . (13) 

By substitution in equation (12) of series in the form of equation (13), 


P(T,AV,p 0 ,r) -por^AY 


* - (^) 


AY 2 + 


AY 3 . . . (14) 
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is obtained. When T is close to zero, is small and only the first 
term of the series in equation (14) need be used. As T increases, the 

terms of higher order must be considered and r^-l becomes more nearly 

equal to rT, because r is greater than 1. Equation (14) may there- 
fore be rewritten as 



The first step in the derivation of equation (10) is the proof 
that Pjj, the mode value of P corresponding to T M for one set 
of values of the parameters p Q and AY, is equal to P M corresponding 
to another set of values . Two sets of values • of Pq and AV are 
indicated by adding the subscripts 1 and 2, and the corresponding 
variables are similarly indicated. With equation (3) in the form 

log (log r) - log p n AY - log 0.4343 

T «* U6) 

with values AV^ and AVg of the parameter AV, with values Pq ^ and 

Pq,2 of the parameter p Q , and with r maintained constant, the equation 

T 1,M t 2,M 

r p 0,l AV 1 - r p 0,2 AV 2 . ( 1? ) 

is obtained by subtraction of the equation for T 2 M from that for T x M . 
The mode probabilities by equation (15) are * 
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.and 



Therefore, by equation (17), the separate terms in the series for P 1 M 
and ,M are equal and 


P - P 

1,M " r 2,M 

Because for a given value of r the mode probability is the 

same for all values of p Q and AV, the probability function P ranges 
from zero to P M for all values of ' p Q and AV. To the value of the 
probability P]_ corresponding to the values p 0 AV^, and T^, there 
are values Pq } 2> AVg, and Tg such that 


From equation (15) series expansions for P-^ and Pg are obtained and 
substituted in this equation so that 

rTl PQ,1 AV 1 " p 0,2 AV 2 = 


(18) 

For values of T close to zero, the probabilities Pj and P 2 are given 

to a good approximation by the first term of the series in equation (15), 
and the right-hand member of equation (18) is therefore negligible. As T 
increases, the terms of second and higher order in equation (15) become 
appreciable, and the individual terms on the right-hand side of equa- 
tion (10) are no longer negligible. But as T approaches T M , the 

bracketed differences on the right-hand side of equation (18) approach 
zero according to equation (17). The differences on the right-hand side 


iik 

r-l 


( 00,1 « 1> 2 - 


(A 

r-l 


(pQ,2 AV 2 y 






2 (P0,1 AVi)2 - ( P 0 ^2 A ^2) 2 
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of equation (18) are therefore negligible in the range of T from 0 
to a value somewhat greater than T^. Therefore, because P becomes 

negligible for T just a little greater than (see fig. 3), the 

equation 


L p AV = r T2 p AV 

0,1 1 0,2 2 


(19) 


is a good approximation in the range of T where P is of appreciable 
magnitude. If the log of equation (19) is taken, this equation may be 
rewritten as a transformation equation, which yields the displacement 
of a probability curve in the T coordinate as 




l0g ^ p 0,l AV l/ p 0,2 AV 2^ 

log r 


( 10 ) 


Calculations of P for various values of p Q and AV from equa- 
tion (2) confirm the approximate validity of equation (10). 
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TABLE I - COMPARISON OF MODES AND AVERAGES OF 
FREQUENCY -DISTRIBUTION CURVES 


Diameter 
of hemi- 
spherical 
drops 
(microns) 

Mode 

(°c) 

Average 

(°c) 

( Mode ) - ( average ) 

(°c) 

^-Experimental values 

805 

-16.6 

-17.4 

0.8 

161 

-23.3 

-24.2 

.9 

69 

-25.5 

-25.2 

-.3 

34.5 

-24.4 

-26.2 

1.8 

8.75 

-32.2 

-32.1 

-.1 

Theoretical values 

1564 

-15.8 

-15.2 

-0.6 

156.4 

-23.0 

-22.6 

-.4 

72.6 

-25.5 

-24.9 

-.6 

33.70 

-27.7 

'-27.3 

-.4 

7.26 

-32.7 

-32.1 

-.6 


iata from reference 2. 


TABLE II - VALUES OF p Q AND CORRESPONDING WEIGHTS 
USED IN MAKING COMPOSITE CURVE 


p o 

Weight 

(cm’3) 


1.0 x io" 3 

2 

2.64 X 10~ 4 

4 

1.05 X 10" 4 

2 

5.7 X 10" 5 

2 

3.7 x 10' 5 

2 

2.6 X 10" 5 

1 







Figure 1. - Probability curves for various values of r computed from exact equation with 
1-AV/V = (n-l)/n = 1 - 10~ 8 , AV = 10 -3 V = 10 5 ; and P n = 2.0X10- 4 . 
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Figure 2. - Straight line fitted to da^a from reference 2 for diameter D against average freezing temperature 
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10 0 -10 -20 
Freezing temperature, °F 


(a) Diameter, D, 805 microns. 

Figure 3* - Comparison of theoretical and experimental probability curves. 
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Freezing temperature, °F 


(b) Diameter, D, 161 microns. 

Figure 3. - Continued* Comparison of' theoretical and experimental 

probability curves. 
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Figure 3. 


1 1 : i 1 

0 -10 -20 -30 

Freezing temperature, °F 

(c) Diameter, D, 69 microns. 

Continued, Comparison of theoretical and. experimental 
probability curves. 
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I 1 1 L_ I I 

10 0 -10 -20 -30 -40 

Freezing temperature, 

(d) Diameter, D, 34.5 microns 

Figure 3. - Continued. Comparison of theoretical and experimental 

probability curves. 
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10 -14 , -18 -22 -26 -30 -34 .-38 

Freezing, temperature , °C 


i : — 1 I 1 ' i l 

10 0 -10 ' -20 -30 -40 

Freezing temperature; °F 

(e) Diameter, D, 8.75 microns. 

Figure 3. - Concluded. Comparison of theoretical and experimental 

probability curves. 



NAC A- Langley - 12-6-50 - 1050 


Figure 4. - Comparison of standard deviations of theoretical distribution curves with those of 






